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Abstract 

Nanocomposites made of polymer films embedding silver nanoparticles were pre¬ 
pared by thermal annealing of poly-(vinyl) alcohol films containing AgN0 3 . Low 
(2.5 % w:w) and high (25 % w:w) doping concentration of silver nitrate were 
considered as well as their effect on the optical properties of thin (30 nm) and 
thick (300 nm and more) films. The topography and the optical properties (re¬ 
fractive index n and extinction coefficient k) of such films were studied by ato¬ 
mic force microscopy and spectroscopic ellipsometry. For a given doping level, 
the parameters of the surface plasmon-polariton resonance (amplitude, position 
and width) were shown to be thickness-dependent. Multivariate statistical analy¬ 
sis techniques (principal component analysis and support vector machines) were 
used to explain the differences in the optical behavior of the thick and thin films. 

Keywords: Spectroscopic ellipsometry, Silver nanoparticles, Polyvinyl alcohol, 
Thermal annealing, Multivariate analysis, Optical properties 


1. Introduction 

Nano-objects and metal nanoparticles (NPs) play a central role in the deve¬ 
lopment of nanotechnology-based optical devices. The optical properties of metal 
NPs originate from the collective oscillations of their conduction electrons. These 
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are termed localized plasmons or surface plasmon polariton resonances (SPPR) 
when the free electrons are excited by the electromagnetic field associated to the 
propagation of light ELS . These optical properties are strongly influenced by the 
shape and size of the NPs but also by the dielectric properties of their environment, 
in particular when the NPs are embedded in a dielectric matrix such as polymers, 
silicon dioxide or titanium dioxide layers [3j]. A wide range of experimental me¬ 
thods is available for the synthesis of such materials. Besides the methods invol¬ 
ving the synthesis of the NPs in a liquid medium, with or without further coating 
to prevent their aggregation (see e.g. fl-0] and the references therein) and their 
dispersion in a solid phase, NPs can be synthesized in situ following the irradia¬ 
tion or the thermal annealing of the solid phase. Although the first approach has an 
evident advantage of leading to the synthesis of NPs with wide ranges of shapes, 
size ..., the synthesis methods of the second category are usually simpler ("one-pot 
synthesis") and avoid a concentration step of the NPs before their dispersion in the 
matrix [@]. In particular, in the case of ’soft’ matrices such as polymers, the repro¬ 
ducibility of the NPs distributions in size, shape ... remains difficult to achieve and 
slight variations in the experimental conditions leads to different results eviden¬ 
ced by modifications of the plasmon resonance band parameters. Moreover, the 
detailed mechanism linking the optical properties of the nanocomposite to their 
structural parameters is today not fully understood although being the subject of 
an increasing number of publications. 

The basis principle of these in situ synthesis methods is the use of one of 
the matrix components as a reducing agent of a metal salt. It is well known that 
polymers are excellent host materials for metal or semiconductor NPs [jfi!] because 
polymers such as poly-(vinyl) pyrrolidone (PVP) or poly-(vinyl) alcohol (PVA) 
can not only play the role of reducing agent but also contribute to the stabilization 
of the NPs. In the latter case, the energy used to promote the reduction of the metal 
cations originates from a thermal annealing process. In particular, the reduction 
step of the silver cations uses electrons originating from the —OH groups of the 
PVA molecules 


110°C 


A g + + PVA Ag —PVA 


60 min 


(i) 


Interest for such Ag/PVA nanocomposites lays in their optical enhancing pro¬ 
perties [10-12], non-linear optical activity f 13M161 and antibacterial or anti-fouling 
properties [17-20Q- Some recent articles have been devoted to the study of the role 
of the polymer matrix [21] or of the chemical nature of the metal salt [2^, 23ll . 

Usually, the optical properties of plasmonic nanocomposites are probed using 
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reflection or transmission spectrophotometry. In our previous publications [24, 


250, we used spectroscopic ellipsometry (SE) to study the optical properties of 


Ag-doped PVA (Ag-PVA) films produced by thermal annealing. As explained he¬ 
reafter, the method has the advantages of being non-destructive and moreover of 
allowing the simultaneous determination of the film thickness and of its complex 
refractive index values. To date, only a reduced number of studies have been car¬ 
ried out using SE as analysis technique [24-28]. Excellent review articles have 
been recently published by Oates, Wormeester and Arwin on spectroscopic ellip¬ 
sometry studies of plasmon resonances at metal-dielectric interfaces of thin films 
and metallic nanostructures [29|, 30]. In |24, 25], we have shown that, in Ag-PVA 


films, increasing the annealing time contributed to the onset of more NPs of smal¬ 
ler size and that the intensity of SPPR was correlated to the density of NPs [25l]. In 
the mid-infrared spectral domain, the presence in the polymer FTIR spectrum of 
additional peaks located at 1134 cm" 1 and 1036 cm" 1 also suggested the presence 
of bonds between Ag and the polymer matrix [24]. 

In this article, we report on the influence the film thickness and on the silver 
doping of the nanocomposite films on the plasmon resonance parameters. More 
precisely, we compare the resonance parameters in thin and thick films at a given 
Ag-doping level. 


INSERT HERE FIG.Q] 


2. Materials and methods 

2.1. Nanocomposite film preparation 


Polymer films were prepared according to the experimental procedure descri¬ 
bed in our previous publications [24, [25] and based on the original method of 
Porel and co-workers 011 32]. A 10% (w:w) poly(vinyl alcohol) solution in Mil- 
liQ water was prepared as a stock solution (PVA, MW = 13,000 - 23,000, 98% 
hydrolyzed, Sigma-Aldrich). Dissolution of the polymer was achieved by heating 
the solution at reflux at 85°C. AgN0 3 (99.99%, Sigma-Aldrich) was added to the 
polymer solution after cooling to room temperature to obtain a silver concentra¬ 
tion of 2.5% or 25%. Percentage corresponds to the w:w silver-to-polymer ratio. 
Required amount of water was added to the polymer stock solution to obtain 2% 
and 8% polymer solutions. The silver-modified polymer solution was spin-coated 
on piranha-cleaned (100) silicon wafers fragments (ACM, France). Silicon sub¬ 
strates were chosen because of their high optical contrast with polymer films and 
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their low intrinsic roughness. Depending on the coating conditions, the final thick¬ 
ness of the dried film ranged from 18-25 nm (2% PVA, 60 sec at 6000 rpm) to 300 
nm and above (8% PVA, 60 sec at 1600 rpm). After being dried in open air atmos¬ 
phere, the coatings were annealed 60 min at 110°C before their optical properties 
being measured by spectroscopic ellipsometry (SE) and their topography imaged 
by atomic force microscopy (AFM). Contrary to the case of glass substrates in 
which, by increasing the temperature, the coatings became yellowish-brown in 
color, in the case of silicon substrates, the color of the samples is mainly domi¬ 
nated by the interferometric color and the contribution of the plasmon resonance 
to the color change can hardly be seen by naked eyes. For comparison, reference 
PVA films were prepared according to the same method but omitting the AgN0 3 
in the polymer solution. 


2.2. Surface characterization 

The topography of the films was studied with a Park XE70 AFM (Park Sys¬ 
tems Corp., Korea), operated in air in intermittent contact mode with commercial 
ACTA tips (resonance frequency: 309 kHz). Areas of 5 /;m x 5 //m, 2 //m x 2 //m 
and 1 /im x 1 //m were imaged for every sample. The resolution of the images was 
256 pixels x 256 pixels. The data were processed with the Gwyddion 2.34 SPM 
images processing software ( http://gwyddion.net /). The roughness of the samples 
was characterized by the average surface roughness parameter (S a ) and by the 
root-mean-square surface roughness parameter (S q ) defined by 
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M and V are length and the width of the image and z hJ is the height of pixel 
(i , j) above the average plane of the sample [33]. 


2.3. Ellipsometry 

Spectroscopic ellipsometry (SE) is an optical analysis technique based on the 
change of polarization of light after reflection of an incident light beam on a sub¬ 
strate or on a stratified sample L34j,|35]. Due to its non-destructive nature, SE is a 
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powerful technique to analyze at once the optical properties of materials and the 
thickness of the different layers of the sample. More specifically, one measures at 
a given wavelength the ellipticity p which is defined by 

p = — = tan \ke* A (4) 

r s 

where r p and r s are the complex reflectance coefficients of the p— and s —components 
of the incident light, respectively, T and A are the ellipsometric angles. They are 
explicitly defined by 


tan 'k = and A = 5 P — S s (5) 

\ r s\ 

where 5 P and 5 S are the phase-shifts undergone by the p— and the s —components 
of the incident wave during the reflection. The angle A therefore measures the re¬ 
lative phase-shift between the perpendicular components of the propagating light 
wave. 

Except in a very limited number of cases, 'k and A cannot be directly conver¬ 
ted into refractive index and/or film thickness. Indeed, the ellipticity is a complex 
quantity which is function of the angle of incidence, of the optical properties of 
the materials and of the thickness of the layers. Taking into account the disper¬ 
sion of the dielectric function, Equation |4] has to be inverted in most cases using 
numerical methods, as suming specular reflection of light at ideal planar interfaces. 

The ellipticity of the samples was measured using a SOPRA GES5E rotating 
polariser spectroscopic ellipsometer with a 1024 channels CCD spectrograph al¬ 
lowing fast and accurate measurements in the range 1.37 to 4.96 eV. The angle of 
incidence (AOI) was set to 70 deg. and the ellipsometer operated in parallel beam 
configuration. SE results were processed using the SOPRA Winelli II software. 

A one-layer Cauchy model was chosen to represent the optical properties of 
the Ag-PVA films in the transparent range and a Lorentzian oscillator was added 
to account for the localized absorption of the plasmon resonance in visible range 

Q : 
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The specific contribution of the plasmon resonance to the wavelength-dependent 
complex dielectric function e(A) is given by 


e(A) — e r (A) + zcj( A) 


(V) 
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where A denotes the wavelength. A 0 is the resonance wavelength of the oscilla¬ 
tor, A its oscillator strength and T 0 its full-width at half maximum (FWHM) and is 
hereafter referenced to as the ’width’ of the resonance, represents the contribu¬ 
tions to e r from the resonances at wavelengths much greater than the measurable 
frequency range. 


2.4. Multivariate statistical analysis 

The multivariate analysis of the data was achieved using the R statistical soft¬ 
ware 1361. The FactoMineR package lf37ll and the el071 package [38] were used 
to carry out the principal component analysis (PCA) of the resonance peaks para¬ 
meters and to carry out the classification tasks, respectively. PCA is a multivariate 
statistical method pioneered by Pearson [39] usually used to reduce the dimen¬ 
sion of a data matrix keeping the maximum variance in the reduced variable set 
[40]. Support vector machines (SVM) are statistical classification algorithms used 
to find the best classifier (i.e. separator) between two sets of (eventually overlap¬ 
ping) data points [41]. Levene test and Kruskal-Wallis post hoc comparisons were 
done using the agricolae package [42]. 


3. Results and discussion 

3.1. Topography of the nanocomposite films 

To improve our investigation at the nanoscale, the topography of the surface 
was measured by AFM in intermittent contact mode after the annealing of the 
film at 110°C during 60 min. Fig.[l]4 and B represent the topography and phase 
images (1 pm x 1 pm) of a typical 30 nm-thick film. The ratio of Ag/PVA was 
25% (w:w). The nanoparticles of silver were clearly observed by a variation of 
the phase value and a local variation of the height. The first observation of these 
images shows that the size of the particles is on average 10 nm. In Fig. [TJC and 
D, we represented the topography and the phase images corresponding to a 300 
nm-thick film at the same Ag/PVA ratio. The size of nanoparticles seems to be 
smaller (on average 5 — 6 nm). 
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S a and S q were calculated with Gwyddion for doped and non-doped represen¬ 
tative samples (Table [[])• The measurements show an significant difference bet¬ 
ween doped and non-doped film in both cases (thick or thin films). As seen from 
both the average and the RMS roughness values, the growth of the nanoparticles 
induced by the annealing of the film controls the roughness of the surface. If the 
thick and the thin samples are compared to each other either in phase or in topo¬ 
graphy, one can see that the particles embedded at the surface of the thick sample 
are smaller than the particles in the thin one. These features are confirmed by exa¬ 
mining the topographic profiles of our films. In Fig. [2l the dashed and the plain 
lines correspond to doped and the undoped samples, respectively. In both case, the 
roughness of the doped surface is higher. 

INSERT HERE TABLE [D 
INSERT HERE FIG.IH 


3.2. Optical properties 

The purpose of this study was to compare the optical properties of thick and 
thin films embedding silver NPs grown in situ from the same polymer solution and 
to determine whether the organization of the NPs was equivalent in both cases. To 
that purpose, the complex refractive index was calculated from SE measurements, 
i.e. from the 'k and A spectra. 

INSERT HERE FIG.0 

One of main advantages of SE over conventional spectrophotometric methods 
is that it allows a simultaneous determination of the optical properties and of the 
thickness of the layers and not only the measurement of absorption or of trans¬ 
mission spectra. Fig. [3] represents the SE data corresponding to a 25.4 nm-thick 
film (Fig. [3jA) and to a 293.9 nm-thick film (Fig. [3^). Instead of representing the 
\k and A ellipsometric angles, we have chosen to present the a = cos(2\k) and 
j3 = sin(2\k) cos(A) Fourier coefficients, which are closer to the experimental 
quantities measured by the ellipsometer. With respect to reference films of equi¬ 
valent thickness, the ellipsometric data present spectra distortion in the 3.0-3.5 
eV range. This is attributed to the localized plasmon resonance associated to the 
silver NPs. This effect is more noticeable for the thick films (Fig. [3j3): one can 
see three peaks in the o-spcctrum (filled symbols) but the second is considerably 
damped due to the energy absorption at the resonance. The data generated by the 
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optimized optical model are also presented in Fig. [3] (dashed lines). The results 
generated by the optical model are in close agreement with the experimental data: 
the residual x 2 values are respectively equal to 2.7 10 -4 and 2.9 10 -3 for the thin 
and thick films. The higher discrepancy between experiments and modeling for 
the thick film is due to a slight thickness inhomogeneity of the samples and to a 
possible distribution of the NPs in the direction perpendicular to the substrate sur¬ 
face that would result is some minor optical anisotropy. This possible cause has 
not been further considered in this study. 

As previously described [24], to model the optical behavior of our samples, 
a one-layer Cauchy model was chosen to represent the optical properties of the 
polymer matrix and a Lorentzian oscillator was added to account for the localized 
absorption of the plasmon resonance in the visible range. The effect of the pre¬ 
sence of the Lorentzian oscillator on the model is equivalent to using the Maxwell- 
Gamett effective medium approximation (MG-EMA) in a view of obtaining the 


dielectric function of the nanocomposite films 11301 . I43L14411 . With this approach, 


two problems are implicitly solved : the dependence of the dielectric function of 
the silver NPs and the generalization of the MG-EMA to take into account the 
size and shape distribution of the NPs. These two problems can easily be solved 


for spherical NPs [lj, 28 


45f] but require more complex approaches for non sphe¬ 


rical particles (see e.g. [2, |46|, 470). Moreover, these advanced models have to be 


supported by complementary experiments based e.g. on the quantitative analysis 
of transmission electron microscopy images to obtain the size distribution of the 
NPs. Also is avoided the influence of the doping level of the films: at high doping 
levels, it is expected that the nanocomposite films cannot be considered anymore 
as diluted solutions of NPs in a polymer matrix. The analysis of the optical res¬ 
ponse of nanocomposite films containing a high concentration of metal nanopar¬ 
ticles is more relevant of the metal island model [48, ;49(] which has been recently 
generalized by Wormeester and co-workers |[50l i5 lh. Adding roughness or void 
inclusions did not significantly improved the quality of the fits, even for the thin 
samples. The use of the Cauchy model for the polymer matrix is justified by the 
monotonous increase of the refractive index of the polymer when going from the 
infrared to the UV, as expected for dielectric materials far from the fundamental 
absorption edge. 


INSERT HERE FIG. [4] 

The optical properties of the NPs-doped films are represented in Fig.[4}'\. For 
each type of film, the refractive index n and the extinction coefficient k of two 
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different samples are represented to show the reproducibility of the analysis and 
indirectly, of the sample preparation. The presence of the silver NPs considerably 
modifies the optical properties in the doped polymer films: the presence of the ab¬ 
sorption peak in the 2.5 - 3.5 eV range of the extinction coefficient spectrum (Fig. 
|4}3) induces a large oscillation in the refractive index value (Fig. |4K) because of 
the Kramers-Kronig consistency of the optical properties and of the Lorentzian 
model of the peak. The perturbation of the refractive index occurs over the 1.5 
- 4.5 eV range although the absorption peak is limited and the refractive index 
values are higher than the ones of a pure polymer film (see e.g. Fig. 5 in refe¬ 
rence [250). The optical properties of thin and thick films are obviously different. 
The optimized parameters of the plasmon resonance peak are given in Tabled In 
thin films, the amplitude of the resonance is higher, the position of the resonance 
smaller and the peaks wider than in thick ones. 

INSERT HERE TABLE |2] 

Data presented in the Supporting information section (Fig.lSTT) show that after 
an annealing time of 40 min. at 110°C the optical properties of the film are stable 
and that all the Ag + ions are reduced. 

3.3. Multivariate analysis 

Once the optical properties of the nanocomposites are determined and remem¬ 
bering that they are normally decorrelated from the thickness of the film, we consi¬ 
dered the dispersion of the data cloud, the optical response of each nanocomposite 
being characterized by 3 parameters (A, A 0 and r 0 ). 

Four categorical variables are defined as a function of the thickness of the film 
(’T’ : thick; ’t’ : thin) and of the doping level in silver (’H’ : high (25%); T : 
low (2.5 %)). In the following, these categorical variables will be referred to as 
the class of the sample. 

3.3.1. Distribution of the resonance parameters 

On the basis of the experimental conditions, the 90 samples considered for 
this study are distributed among 4 classes relatively equivalent sizes. Each of the 
four categorical variables (i.e. nanocomposite classes) are represented by a dif¬ 
ferent color. Let us first consider the link between the thickness of the film and the 
strength of the plasmon band A (Fig. [5]). One can easily see that the data cloud is 
subdivided into four distinct classes, each of them being assigned to one type of 
nanocomposite. 


9 


INSERT HERE FIG.© 
INSERT HERE FIG.© 


As one can see from the scatter plot matrix of the resonance parameters (Fig. 
[6]), the subclouds are much more overlapping, particularly in the case of the T 0 - 
A 0 plot. In the first row of the scatterplot matrix, different behaviors can be de¬ 
tected between the nanocomposites prepared from solution with a high AgN0 3 
concentration (red and green circles) and those prepared from diluted AgN0 3 so¬ 
lution (blue and yellow circles). 

INSERT HERE TABLE [3] 

This global picture of the data cloud is reinforced by the statistical dispersion 
of the resonance parameters in terms of box-and-whiskers plots (Fig. [7]). Heteros- 
cedasticity of the data was checked using a robust version of Levene test. The re¬ 
sults show the null hypothesis H 0 (equality of all the variances) has to be rejected 
at a confidence level better than 1% (data not shown) for all the resonance para¬ 
meters. Although classical analysis of variance (ANOVA) methods are relatively 
robust with respect to the non-normal distribution of the variables, differences 
in variances are more complicated to handle. As a rule of the thumb, classical 
ANOVA can be used if the ratio of the largest to the smallest variance is < 4. This 
case is not met by our data. Non parametric tests are therefore required. Results 
issued from the Kruskal-Wallis test are presented in Table |4j In a given row of the 
table, classes represented by the same letter do not statistically differ form each 
other. 


INSERT HERE FIG. ID 

As expected, the strength of the oscillator A is strongly dependent on the silver 
content of the film (Fig. |7]A). Two trends emerge from the peak position disper¬ 
sion : increasing the thickness of the films (t —» T) slightly blue-shifts the reso¬ 
nance wavelength, as well as increasing the amount of Ag + present in the films 
(l —> H ) (Fig. |7}3). The evolution of the resonance width is more complicated 
(Fig. |7p). For that reason, the behaviors of the films with low and high silver 
contents have to be separately described in the remaining parts of this paper. 

INSERT HERE TABLES) 
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3.3.2. Correlation between the plasmon resonance parameters 

To evaluate the correlation between the plasmon resonance parameters ( A , 
L 0 and Gamma) and their link with the experimental variables (thickness of the 
film and silver content), we performed a principal components analysis to find the 
principal axes of the data cloud. Thin and thick films were separately considered. 
Taking into account the reduced number of active variables (3), the percentage of 
the variance expressed in the first factorial plan is high (Table [5]). Additional va¬ 
riables (thick : film thickness and SilverC : silver concentration) are also projected 
in the first factorial plan and are represented in blue in Fig. [8] For both the thin 
and the thick films, the silver concentration is, as expected, strongly correlated to 
the strength of the oscillator A. More surprising is the fact that the thickness is 
not correlated with the plasmon resonance parameters for thin films (Fig. [8[ top) 
and appears clearly to be anti-correlated with the width of the resonance for the 
thicker ones (Fig. [8] bottom). 

INSERT HERE TABLE 0 

INSERT HERE FIG.! 

3.4. Classification in the r 0 — A 0 plane 

As we have seen from the PC analysis, thick and thin films differently behave 
when we consider the link between the resonance width (T 0 ) and the position of 
the plasmonic band (A 0 ). We therefore attempted to classify the films as a function 
of the silver content of the coated solution (2.5% or 25%). The classification task 
was performed for the thin and the thick films separately (Fig. [9]). 

When the data are presented in the r 0 — A 0 plane, one can see that the set 
of data points corresponding to the highly doped films (circles) has only a slight 
overlap with the cloud of the films prepared from the 2.5% solution (triangles). 
To perform the classification task in a quantitative way, we used a SVM algorithm 
[38] derived from the original classification algorithm proposed by Vapnik fiH. 

SVMs are supervised learning algorithms used to find the best classifier (i.e. 
separator) between two (or more) data sets, overlapping or not. In the simplest 
cases (two non-overlapping data clouds in a 2-dimensions space), the algorithm 
determines the line whose position maximizes the distance to the nearest data 
points belonging to each class (i.e. the margin). For that reason, SVMs are also 
named large margin classifiers. In more complicated cases, non-linear classifiers 
can be used. One important characteristics of the SVMs has to be kept in mind : the 
equation of the separator can only be explicitly determined in the case of a linear 
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separator. In other cases, the equation of the separator is not known because the 
classification task is performed is a derived space (usually of a higher dimension 
than the data space) on the basis of an inner product only Eh. Besides the choice 
of the kernel (linear or not), the classifier has to be tuned with respect to the penalty 
applied when an event is misclassified. In each case, the classifier has been tuned 
using a 10 x cross-validation process and varying the penalty C in the range 0.01 
to 10. Best results were obtained with a cost C for the misclassification equal to 
0.15 and 0.40 for the thin and thick films respectively. Optimum results for the 
classifiers are represented in Fig. [9] The figure shows that, when applied to our 
data, SVMs allowed us to discriminate between the high- and low-doping level 
films in a very efficient way : only 3 films over 33 were misclassified in the thin 
films case and 3 over 54 in the thick films one. The most striking feature of the 
classification is the behavior of the classifiers : when comparing panels (A) and 
(B) in Fig.lH one can easily see that their slope have opposite signs. The classifiers 
slopes are respectively -0.36 (thick films) and +0.30 (thin films). 

INSERT HERE FIG. [9] 

These different optical behaviors are probably strongly related to the growth 
process of the NPs in films of different thicknesses. In thin films, NPs grow in 
a 2D-like matrix while in thick films the NPs rapidly form a 3D structure. In 
forthcoming papers, we will investigate using AFM the real-time growth of the 
silver NPs in thin and thick films as well as their final spacial distribution based 
on electron microscopy results. 

4. Conclusions 

In summary, we have studied the optical properties of nanocomposite contai¬ 
ning in situ grown nanoparticles. In particular, the study has been focused on the 
influence of the thickness of the film and of the silver concentration in the po¬ 
lymer matrix on the optical properties of the nanocomposite. We experimentally 
observed a different behavior of the thick and thin films as a function of the sil¬ 
ver concentration. Starting from spectroscopic ellipsometry results, the statistical 
distribution of the plasmon resonance parameters has been analyzed using princi¬ 
pal component analysis and large margin classifiers. Both techniques lead to the 
conclusion that at a given Ag + doping level, thin and thick films behave differently, 
probably due to a different spatial distribution of the NPs in the polymer matrix. 
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FIGURE 1: Topography (left) and phase (right) AFM images in intermittent contact mode of Ag- 
PVA film doped with 25% AgN0 3 (w:w). A,B : ~ 30 nm-thick film; C,D : ~ 300 nm-thick film. 
Image size: 1 /xm x 1 /xm (256 x 256 pixels). 
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FIGURE 3: Ellipsometric spectra of Ag nanoparticles-doped PVA films (Ag/PVA ratio: 25 % w:w): 
A, thin films (thickness: 25.4 nm); B, thick film (thickness: 293.9 nm). Experimental data: a = 
cos(2'I') (filled circles) and /3 = sin(2'I') cos(A) (open circles). Dashed lines: optimized results 
from the optical model. 
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Figure 4: Optical properties of thin (plain lines) and thick (dashed lines) silver NPs-doped PVA 
films (Ag/PVA ratio: 25% w:w): A, refractive index n; B, extinction coefficient k. 
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Figure 5: Scatterplot of the thickness of the film (in /im) versus the strength of the oscillator. 
Color coding of the categorical variables are given in the legend. 
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FIGURE 6: Scatterplot matrix of the resonance parameters. Colour coding of the categorical va¬ 
riables are given in the legend. Units for peak position ("LO" or Ao) and peak width ("Gamma" or 
Tq) are micrometers. 
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FIGURE 7: Box-and-whiskers plots of the resonance peaks parameters : peak height (A), peak 
position (L0) and peak width (Gamma). Categorical class definitions are given in the text. 
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FIGURE 8: Correlation circles in the first factorial plane : (top) thin films and (bottom) thick films. 
Black : Variables used to perfrom the anaysis; Blue : additional projected variables (thickness and 
siver content of the coated solution). 
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Figure 9: Classification of the weakly and highly Ag-doped film in the To - Ao plane. (A) Thick 
films (B) Thin films (Ag-PVA ratio: 25% (circles) and 2.5% (triangles); Open symbols: support 
vectors; Lines : optimum classifier (plain) and margin (dashed). 
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Thick film 

Thin film 



Doped Undoped 

Doped Undoped 

5a 

(nm) 

0.66 

0.40 

0.59 

0.32 

s q 

(nm) 

0.85 

0.52 

0.77 

0.42 


Table 1: Roughness parameters of the doped (Ag/PVA = 25%) and undoped (pure PVA) films 


Sample 

d (nm) 

A 

A 0 (nm) 

To (nm) 

Thin 

23.4 ± 0.2 

0.145 ±0.6 

414.2 ±0.7 

67.6 ±2.9 


25.4 ± 0.3 

0.133 ±0.5 

415.6 ±0.6 

69.0 ± 2.6 

Thick 

305.9 ± 1.7 

0.117 ± 0.2 

405.4 ± 0.7 

47.3 ± 1.6 


293.4 ± 1.7 

0.118 ±0.2 

409.5 ± 0.6 

49.2 ± 1.5 


Table 2: Typical parameters of the plasmon absorption peak (A: amplitude of the absorption 
peak; Ao: position of the resonance; To: width of the resonance) as a function of the film thickness 
for highly doped PVA films (Ag/PVA ratio: 25% w:w). Data correspond to the optical properties 
presented in Fig.Q] 


class 

N 

thick 

Means 

A L0 

Gamma 

thick 

Standard deviations 

A L0 Gamma 

tH 

17 

0.028 

0.125 

0.421 

0.095 

0.014 

0.050 

0.006 

0.046 

TH 

27 

0.311 

0.102 

0.410 

0.064 

0.082 

0.018 

0.005 

0.018 

tl 

16 

0.030 

0.013 

0.437 

0.064 

0.019 

0.014 

0.014 

0.036 

T1 

30 

0.254 

0.028 

0.419 

0.127 

0.076 

0.024 

0.013 

0.056 


Table 3; Class sizes ( N ), means and standard deviations of the resonance parameters as a function 
of the nanocomposite type. Units for thickness ("thick"), peak position ("LO" or Ao) and peak width 
("Gamma" or To) are micrometers. 


Parameter 

tH 

Class 

TH tl 

Tl 

A 

a 

a 

b 

c 

Ao 

b 

c 

a 

b 

r 0 

b 

c 

c 

a 


Table 4: Summary of the Kruskal-Wallis statistical test for the resonance parameters. In a gi¬ 
ven row, classes represented by the same letter do not statistically differ form each other. To be 
compared with the boxplots in Fig. [7] 
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Samples 

Components 

Eigenvalue 

Variance (%) 

Cum. variance (%) 

Thick 

1 

1.82 

60.83 

60.83 


2 

0.89 

29.62 

90.44 


3 

0.29 

9.56 

100.00 

Thin 

1 

1.72 

57.32 

57.32 


2 

0.85 

28.32 

85.64 


3 

0.43 

14.36 

100.00 


Table 5: Eigenvalues, percentage and cumulative percentage of the variance for the thick and thin 
films 
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Supporting information 

Figure [ST| represents the time evolution of the height of the resonance peak 
observed for Ag-PVA films on glass substrates during annealing at 110°C. The 
films were prepared in the same experimental conditions as those used for the 
’TFT samples on silicon. The absorbance was measured using a Genesys 20 UV- 
Vis spectrometer (Thermo, USA). The spectra (data not shown) were recorded 
from 350 nm to 850 nm. As shown on the figure, the absorbance stabilizes after 
40 min., indicating the complete reduction of the Ag + ions. The fluctuations are 
due to the changes in thickness from one sample to the other. 
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FIGURE S 1: Time evolution of the absorbance at peak maximum ( Annealing temp.: 110°C - Data 
point: median. Error bars: interquartile range, dashed lines: weighed regression lines). 
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